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Abstract 

Ongoing and future wide-field galaxy surveys can be used to locate a number of clusters of 
galaxies with cosmic shear measurement alone. We study constraints on cosmological models 
using statistics of weak lensing selected galaxy clusters. We extend our previous theoretical 
framework to model the statistical properties of clusters in variants of cosmological models as 
well as in the standard ACDM model. Weak lensing selection of clusters does not rely on the 
conventional assumption such as the relation between luminosity and mass and/or hydrostatic 
equilibrium, but a number of observational effects compromise robust identification. We use 
a large set of realistic mock weak-lensing catalogs as well as analytic models to perform a 
Fisher analysis and make forecast for constraining two competing cosmological models, wCDM 
model and f{R) model proposed by Hu & Sawicki, with our lensing statistics. We show that 
weak lensing selected clusters are excellent probe of cosmology when combined with cosmic 
shear power spectrum even in presence of galaxy shape noise and masked regions. With the 
information of weak lensing selected clusters, the precision of cosmological parameter estimate 
can be improved by a factor of 1.6 and -- 8 for wCDM model and f{R) model, respectively. 
Hyper Suprime-Cam survey with sky coverage of 1250 squared degrees can constrain the 
equation of state of dark energy wo with a level of Awq ~0.1. It can also constrain the additional 
scalar degree of freedom in f{R) model with a level of l/^ol 5 x 10“®, when constraints from 
cosmic microwave background measurements are incorporated. Future weak lensing surveys 
with sky coverage of 20,000 squared degrees will place tighter constraints on wq and I/rq] even 
without cosmic microwave background measurements. 


1 INTRODUCTION 

Statistical analyses of large astronomical data set have estab¬ 
lished the standard cosmological model called ACDM model. 
Various statistical methods have been proposed, and have been 
actually applied to real observational data. Well-known ex¬ 
amples include luminosity distance to Type la supernovae 


(Betoule et al. 2014), baryon acoustic oscillation measure¬ 
ment with galaxy clustering (Beutler et al. 2011; Blake et 
al. 2011; Anderson et al. 2014), the anisotropies of cosmic 
microwave background (CMB) (Hinshaw et al. 2013; Planck 
Collaboration et al. 2014), measure of gravitational growth 
by large-scale structure (Sanchez et al. 2012; Beutler et al. 
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2014), and weak gravitational leasing (Kilbinger et al. 2013). 
Although the standard ACDM model is consistent with the re¬ 
sults from a broad range of cosmological analyses, several alter¬ 
native models are still allowed at present. In particular, the ap¬ 
parent cosmic acceleration can be explained not only by ACDM 
model but also by a different class of cosmological models that 
invoke evolving dark energy or modified gravity. Such alter¬ 
native models are basically constructed by the modification of 
Einstein equation: 

G^^ = 8nGT^^. ( 1 ) 

Dark energy models suppose an exotic form of energy in the 
right hand side of Eq. (1), while modified gravity models 
changes the left hand side of Eq. (1) without assuming an un¬ 
known energy. In order to test the two scenarios, it is essential to 
measure the gravitational growth of matter density fluctuations, 
because the modification of gravity could induce characteristic 
clustering patterns to the matter density distribution in the uni¬ 
verse. 

Gravitational leasing is one of the most promising tools to 
probe matter density distribution in the universe. The fore¬ 
ground gravitational field causes small image distortion of dis¬ 
tant galaxies. The small distortion contains, collectively, rich 
cosmological information about the matter distribution. Since 
image distortion induced by gravitational leasing is very small 
in general, we require the statistical analysis in order to extract 
cosmological information from gravitational leasing. Among 
various statistical methods of gravitational leasing, mapping 
matter distribution on a continuous sky is a key, basic method 
(Massey et al. 2007; Van Waerbeke et al. 2013; Chang et 
al. 2015). The reconstruction is purely based on deflection of 
light from sources and thus free of conventional astrophysical 
assumptions. Non-Gaussian features in the reconstructed den¬ 
sity field are produced by non-linear gravitational growth and 
can not be extracted by means of the conventional statistics of 
cosmic shear such as two-point correlation function or power 
spectrum. 

Clusters of galaxies serve as a promising probe of cosmol¬ 
ogy. The number count of clusters is expected to be highly sen¬ 
sitive to growth of matter density perturbations (Lilje 1992), 
whereas correlation analysis of the position of clusters and cos¬ 
mic shear can constrain both the matter density profile and clus¬ 
tering of clusters (Oguri et al. 2012; Okabe et al. 2013;Covone 
et al. 2014). The advantage of weak lensing among various 
techniques is that it does not rely on physical state of the bary- 
onic component in clusters. The finding algorithm with cos¬ 
mic shear analysis is based on reconstruction of matter density 
distribution along a line of sight (Hamana, Takada & Yoshida 
2004; Hennawi & Spergel 2005; Maturi et al. 2005; Marian et 
al. 2012). Reconstructed mass density map is utilized to iden¬ 
tify high density regions in the universe that correspond to mas¬ 


sive collapsed objects such as clusters of galaxies (Miyazaki et 
al. 2007; Schirmer et al. 2007; Shan et al. 2012). 

Previous studies investigated cosmological information in 
number counts of weak lensing selected clusters (Maturi et al. 
2010; Kratochvil, Haiman & May 2010 ; Dietrich & Hartlap 
2010; Yang et al. 2011; Hilbert et al. 2012). In our earlier 
work (Shirasaki et al. 2015, hereafter Paper I), we examined 
other statistics beyond the abundance of weak lensing selected 
clusters. There, we study the statistical property of weak lens¬ 
ing selected clusters using realistic mock weak lensing catalogs 
that incorporate masked regions and shape noise contaminant. 
We developed a theoretical framework based on halo model that 
provide robust predictions for the statistical properties of comic 
shear field and weak lensing selected clusters. In this paper, we 
derive constrain on variants of cosmological models with the 
lensing statistics. We extend our theoretical model to predict 
the lensing statistics for two competing cosmological models 
(dark energy model and modified gravity model). In order to 
realize the realistic situation in galaxy imaging surveys, we use 
two hundreds mock weak lensing catalogs with the proposed 
sky coverage of ongoing Subaru Hyper Suprime-Cam (HSC) 
survey* (see Paper I for details). We derive accurate covariance 
matrices between the lensing statistics using our large mock cat¬ 
alogs. We then make realistic forecast for cosmological con¬ 
straints on two competing scenarios with weak lensing selected 
clusters. 

The paper is organized as follows. In Section 2, we briefly 
describe the basics of two cosmological models called dark en¬ 
ergy model and modified gravity model. There, we describe in 
detail the evolution of background and linear density perturba¬ 
tions in each model. We also summarize the statistical property 
of weak lensing selected clusters and the theoretical model of 
weak lensing statistics in Section 3. By using our theoretical 
model and a large set of mock catalogs, we make forecast for 
constraining on cosmological models by our method in the up¬ 
coming HSC survey in Section 4. Conclusions and discussions 
are summarized in Section 5. 

2 COSMOLOGICAL MODEL 

There exist various extensions to the standard ACDM model. 
We consider two competing models: tnCDM model and f{R) 
model. The former corresponds to a cosmological model based 
on General Relativity with dark energy and the latter represents 
a model with modified gravity. Either model can explain the 
observed acceleration of cosmic expansion al z with appro¬ 
priate parameters. Throughout this paper, we assume a spatially 
flat universe. 

^ http://www.naoj .org/Projects/HSC/index.html 
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2.1 wCDM model 


The expansion of the universe is described by the scale factor 
a{t). We adopt the usual normalization a = 1 at present. In 
General Relativity under a Friedmann Robertson Walker (FRW) 
metric, one can derive the time evolution of a{t) as 


1 da 

Hia) = —— 
^ ' a At 


— FfosGmoa + Gde 


X exp 


—3^ do'(l + WDE(a'))/a' 


1/2 


( 2 ) 


where H{a) is known as the Hubble parameter and Hq — 
lOO/i km Mpc“^ is the present value of Hia). The 
equation of state of dark energy is specified as WDE(a) = 
F’DE(a)/pDE(a) where Pde and poE denote to pressure and 
density of dark energy, respectively. 

In a wCDM model, the linear growth rate of matter density 
fluctuations after matter domination is given by the solution of 
the following equation (Wang & Steinhardt 1998): 


A^g+ [5 
dlna2 ^ [2 


-iUDE(a)HDE(a) 


dff+ 

dlna 



uiDE(a)] HDE(a)5i+ = 0, 


(3) 


where g+(a) = D(a)/a and D{a) represents the linear growth 
rate of matter density perturbations. We can obtain D{a) 
by solving Eq. (3) with the boundary conditions g+{a) — 1 
and d(/+/dlna = 0 at a <C 1. The linear density perturba¬ 
tions are commonly characterized by the linear power spectrum 
P^{k,a). In matter domination, the power spectrum of primor¬ 
dial curvature perturbation is related to that of the matter density 
through Poisson equation as follows: 


47rfc^P^(fc,a) ,, 2 /, ^ 

(2^)3 =^^(^ 0 ) 


r^mo 


D^{a)T^{k) 


ii) 


— l+TTa 


a2.,x 47vk Pn{k) 
- (2^)3 ’ 


(4) 

(5) 


where T{k) is the transfer function of matter density fluctua¬ 
tions and PTi{k) represents power spectrum of curvature per¬ 
turbation. In this paper, we normalize P^ with the value of 
As = A^{ko) at ko = 0.002Mpc“^ as follows in Hinshaw et 
al. (2013). 

Throughout this paper, we consider the simplest model of 
dark energy with constant value of wo- The non-linear gravita¬ 
tional growth of matter density fluctuations in the model have 
been investigated by previous numerical studies. For various 
Wo, a large set of cosmological simulations have been used to 
derive accurate non-linear power spectrum (Takahashi et al. 
2012; Heitmann et al. 2014). Also, Bhattacharya et al. (2011) 
have studied the abundance of dark matter halos for a wide 


range of cosmological parameters in uiCDM model. In this 
work, we pay special attention in order to derive observational 
constraints on wo from non-linear cosmological information. 


2.2 f(R) model 

In f(R) model, the Einstein-Hilbert action is modified by a gen¬ 
eral function of the scalar curvature R, 

'R + f(R) 


Sg= a 


( 6 ) 


IGttG 

The action with Eq. (6) leads the modified Einstein equation as 


+ fRRp.1, “ ( 2 ~ 


= StvGTu 


(7) 


where Jr = d//di?. Assuming a FRW metric, one can deter¬ 
mine the time evolution of the Hubble parameter in f(R) model 
as follows: 


H'^ 



AH 

dlna 



+ ^+H^fRR^ 
SttG 


( 8 ) 


One can also consider the evolution of matter density pertur¬ 
bations in f{R) model. For sub-horizon modes {k >, aH) in the 
quasi-static limit^, the linear growth of matter density perturba¬ 
tions is determined by (Bean et al. 2007) 


dV (3 1 AH\ 

do^ V a H Aa ) 


dfl+ 

do 


X 


3Hmoa ^ 

(Ff/Ho)^(l + /«) 


1- 2Q \ 

2- 3q) 



(9) 


where Hmo is the effective matter density at present time. We 
can specify this effective density Hmo as 


77/(ij) = 77o s Hmoa + Hde 


X exp 


—3/ da'(l+ «)DE(a'))/o' 


1/2 


( 10 ) 


where is given by Eq. (8). The function Q in Eq. (9) is 

given by 

Q{k,a) = -2(-y (11) 

\aj 1 + jR 

Note that the function of Q induces non-trivial scale depen¬ 
dence of the linear growth rate g+{k,a) — D{k,a)/a in f{R) 
model, while the linear growth rate is solely a function of a in 
General Relativity. 

In this paper, we will consider the representative example 
of f{R) model proposed by Hu & Sawicki (2007) (hereafter 
denoted as HS model). 


^ de La Cruz-Dombriz, Dobado & Maroto (2008) have shown that the quasi¬ 
static approximation becomes quite reasonabie for modeis with |/k| <g 1 
today. 
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f{R) = -2A 


7?" 

’ 


( 12 ) 


where A, p and n are free parameters in the model. For R^ , 
one can approximate the function of f{R) as follows: 


f{R)^-2A- 


fm 

n R" ’ 


(13) 


where Ro is defined by the present scalar curvature of the 
background space-time and fno = —2Ap^/Ro — fR{Ro). In 
the HS model with Ifm] <C 1, the background expansion be¬ 
haves similarly to the standard ACDM model. In practice, for 
I /ho I 10 “^, geometric tests such as measurement of super¬ 
novas cannot distinguish with HS model and the ACDM model 
(Martinelli et al. 2012). Nevertheless, measurement of gravi¬ 
tational growth would be helpful to constrain on HS model if 
the scale dependence of growth rate as shown in Eq. (9) is de¬ 
tected. The non-linear gravitational growth in HS model have 
been studied with cosmological simulations (Oyaizu, Lima & 
Hu 2008; Schmidt et al. 2009a; He, Li & Jing 2013; Zhao 
2014). These studies suggest that statistics of galaxy groups 
or clusters provide meaningful information about the modifica¬ 
tion of gravity. Therefore, a combination of statistics of cosmic 
shear and galaxy clusters can offer an interesting probe of HS 
model. In the following, we focus on the case of n = 1 for which 
many numerical studies have been performed. 


3 Weak lensing 

Here, we summarize basics of weak gravitational lensing and 
the finding algorithm of galaxy clusters with weak lensing mea¬ 
surement. Further details of the statistical properties of weak 
lensing selected clusters are found in Paper 1. 


3.1 Basics 

When considering the observed position of a source object as 
6 and the true position as (3, one can express the distortion of 
image of a source object by the following 2D matrix: 

(14) 

where k is convergence and 7 is shear. 

Let us consider the case of General Relativity. In this case, 
one can relate each component of Aij to the second derivative of 
the gravitational potential as follows (Bartelmann & Schneider 
2001; Munshi et al. 2008); 


_ f 1 - K - 71 -72 

dOf —72 1 — K + 71 


Aij — Sij "Fij, 


'^‘7 = 4 / ^xg{x,x)didj^{x), 

^ Jo 


g{x,x) = 


r{x-x')r{x') 

r(x) 


(15) 

(16) 

(17) 


where x is the comoving distance and r{x) represents the co¬ 


moving angular diameter distance. Gravitational potential <F 
can be related to matter density perturbation 5 according to 
Poisson equation. Therefore, convergence can be expressed as 
the weighted integral of <5 along the line of sight; 

«=|(^) Xl^o ( ( 18 ) 

J 0 

Next, we consider weak gravitational lensing effect in HS 
model. When adopting the Newtonian gauge, one can express 
the line element as 


ds^ = a{rtf [- (1 -f 2$) Arf -f (1 - 2^) Ax^] , (19) 


where 77 is the conformal time. In this metric, the deflection 
angle by gravitational lensing can be written as (Bartelmann & 
Schneider 2001; Munshi et al. 2008), 



where corresponds to the perpendicular component of the 
gradient along a line of sight for proper distance and A£ repre¬ 
sents the line integral in terms of proper distance. Considering 
the modified Einsfein equation in HS model with |4>| <C 1, 
14^1 -C 1 , and |/ho| 1 , one can derive the following equations 

of 4> and tp (the derivation is found in e.g., Arnold, Puchwein 
& Springel (2014)); 


V 4>=—— 
a-‘ 3 6 


( 21 ) 

( 22 ) 


where V is the gradient with respect to comoving distance and 
5p and 5R represent the density perturbation and the fluctua¬ 
tion of the scalar curvature R, respectively. Hence, the lens¬ 
ing potential (<I> -h 4')/2 follows the same equation in General 
Relativity as 

=47rG5p. (23) 

These results show that Eqs. (15), (16) and (18) are available in 
the HS model with |/ho| -C 1. 


3.2 Cluster selection 

Weak lensing is a powerful tool to reconstruct the projected mat¬ 
ter density field. The conventional technique for reconstruction 
is based on the smoothed map of cosmic shear. Let us first de¬ 
fine the smoothed convergence field as 


/C(0) = / d>K(6»-(/))[/(</)), 


(24) 


where U is the filter function to be specified below. We adopt 
the compensated Gaussian Alter for U as 




1 

tv6o 


1 — exp ( — 


el 


(25) 


where 9o represents the boundary of the filter and we set U to 
be zero for 6 > 9o. 
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For a given function of U, the power spectrum of a noise 
convergence field M is given by (Van Waerbeke 2000) 

Pu{l) = ^\U{lt (26) 

^^gal 

where is the rms of the intrinsic source ellipticities rigai rep¬ 
resents the number density of source galaxies, and U is the 
Fourier transform of U. From Eq. (26), we define the moment 
of Af as 

/ r \ 

anois.,i^i Pm Wj ■ (27) 

Throughout this paper, we set ay = 0.4, rigai = 10arcmin“^ and 
assume the source redshift of ^source = 1. These are typical val¬ 
ues for ground-based galaxy imaging surveys (Fteymans et al. 
2012 ). Also, we adopt the smoothing scale of 6a = 5/V81n2 = 
2.12 arcmin and Oo = 30 arcmin. This leads to anoise.o — 0.017. 
The set up is examined in detail in Paper I with numerical sim¬ 
ulations. 

On a smoothed lensing map, convergence peaks with high 
signal-to-noise ratio u = /C/o-noise.o are likely caused by galaxy 
clusters (Hamana, Takada & Yoshida 2004). We thus locate 
high-;^ peaks on a /C map and associate each of them with an 
isolated massive halo along the line of sight. For a dark matter 
halo, we assume the universal NFW density profile (Navarro, 
Frenk & White 1997). We adopt the functional form of the 
concentration parameter in Duffy et al. (2008), 

/ \ - 0.081 

Note that the NFW density profile provides a reasonable flf also 
for halos in HS model with |/ijo| <S 1 (Schmidt et al. 2009a; 
Lombriser et al. 2012) and lhat Ihe corresponding convergence 
Kh can be calculated analytically (Hamana, Takada & Yoshida 
2004). 

In order to predict peak heights in a /C map, we adopt the 
simple assumption that each peak position is exactly at the halo 
center. Under this assumption, the peak height in absence of 
shape noise is given by 

/Cpeak,/i = a J (i^(l)U{(j>;9G,6o)K,h{(t>) +P, (29) 

where a = 0.9 and /? = 0 are found to be in good agree¬ 
ment with numerical simulations as shown in Paper I. The 
actual peak height on a noisy K. map is determined not by 
Eq. (29), but by a probability distribution function (Fan, Shan 
& Liu 2010). The probability distribution function of the mea¬ 
sured peak height /Cpeak.obs with a given /Cpeak,/i is denoted by 
Prob(/Cpeak.obs|Aipeak,/i) in tMs paper. The detailed functional 
form of Prob(/Cpeak.obs|A^peak,/i) IS fouud in Paper I. In Paper 
I, we show that our model of Prob(/Cpeak,obs|A^peak,/i) is quite 
reasonable for peaks with high signal-to-noise ratio in the case 
of 9a ^ 2 arcmin, Oy = 0.4, and ngai 10arcmin“^. 


3.3 Statistics 

We here summarize a set of statistics obtained from weak lens¬ 
ing measurement. Table 1 shows the summary of our model of 
each lensing statistic. 

3.3.1 Convergence power spectrum 

First, we consider convergence power spectrum, which is a 
direct probe of the underlying matter density held. Under 
the Limber approximation^ (Limber 1954; Kaiser 1992 ) and 
Eq. (18), one can calculate the convergence power spectrum as 

where Ps {k) is the three dimensional matter power spectrum, 
Xs is comoving distance to source galaxies and Wn{x} is the 
lensing weight function defined as 

W4X) = I + (31) 

The non-linear gravitational growth of Ps(k) significantly 
affects the amplitude of convergence power spectrum at the 
angular scales less than 1 degree (Jain, Seljak & White 
2000; Hilbert et al. 2009; Sato et al. 2009). Therefore, accu¬ 
rate theoretical prediction of the non-linear matter power spec¬ 
trum is essential for cosmological constraints from weak lens¬ 
ing power spectrum. In order to predict the non-linear evolution 
of Ps{k) for the standard ACDM universe, numerical approach 
based on W-body simulations has been employed extensively 
over the past few decades (Peacock & Dodds 1996; Smith et al. 
2003; Heitmann et al. 2010; Takahashi et al. 2012). In particu¬ 
lar, Takahashi et al. (2012) provides an accurate fitting formula 
of non-linear Ps{k) for various wCDM models. We adopt their 
model in the following. For the HS model of modified grav¬ 
ity, Oyaizu, Lima & Hu (2008); He, Li & ling (2013); Zhao 
(2014) investigated the non-linear evolution of Ps{k). Recently, 
Zhao (2014) derived a new fitting formula of non-linear Ps{k) 
with a large set of numerical simulations. We use the model in 
Zhao (2014) for the HS model. 

3.3.2 Convergence peak count 

We identify local maxima in a smoothed lensing map and match 
each peak with a massive dark matter halo along the line of 
sight. The simple peak count is useful to extract the information 
of the abundance of massive clusters. One can select the lensing 
peaks by its peak height. We define fhe peak signal-to-noise 
ratio by i/ = /Cpeak.obs/cTnoise.o. For a given threshold r-thre, one 
can predict the surface number density of peaks with u > j/thre 
as follows (Hamana, Takada & Yoshida 2004): 

^ The validity of Umber approximation have been discussed in e.g., Jeong, 
Komatsu & Jain (2009). The typical accuracy of Limber approximaion is of 
a level of .C 1 % for ^ > lo. 
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Statistic 

Definition 

Integrand in halo model 

Reference 

Convergence power spectrum 

Eq. (30) 

Matter power spectrum P? 

Takahashi et al. (2012) (wCDM) 
Zhao (2014) (/(P)) 

Convergence peak count iVpeak 

Eq. (32) 

Halo mass function dn/dM 

Bhattacharya et al. (2011) (wCDM) 
Li&Hu (2011) (/(P)) 

Convergence peak auto spectrum Ppp 

Eq. (34) 

Halo mass function dn/dM 
Linear halo bias bh 

Bhattacharya et al. (2011) (wCDM) 
Li&Hu (2011) (/(P)) 

Convergence peak cross spectrum Pp„ 

Eqs. (37) and (38) 

Halo mass function dn/dM 
Linear halo bias bh 

Bhattacharya et al. (2011) (wCDM) 
Li&Hu (2011) (/(P)) 


Table 1. We summarize the elements in our model to derive the lensing statistics presented in this paper. Each column shows, the 
statistical quantity of interest, the definition, the integrand in our model and the reference of fitting formula to compute the integrand. 


Wpeak(r^th re) 


J dzdM 


dV dw 

dzdrt dM 


S{z,M\uthre), 


(32) 


where dn/dM represents the mass function of dark matter halo 
and the volume element is expressed as d^V/dzdfJ = /H{z) 

for a spatially flat universe. In Eq. (32), S{z,M\vthTe) expresses 
the selection function of weak lensing selected clusters for a 
given threshold of uthie- It is given by 


S{z,M\Uthre) = 


L 


d/Cpeak,obs 


^thre^noise,0 


X Prob(/Cpeak,obs | Aipeak,h • 


(33) 


Throughout this paper, we adopt the model of halo mass 
function in Bhattacharya et al. (2011) for wCDM models and 
the model in Li & Hu (2011) for HS models. 


3.3.3 Convergence peak auto spectrum and cross spectrum 

Next, we consider the auto-correlation function of peaks and 
the cross correlation of peaks and convergence. We expect that 
these statistics contain information of dark matter density pro¬ 
file around clusters as well as of the clustering of clusters. 

In Paper I, by using halo model approach, we derive the auto 
power spectrum of weak lensing selected clusters for a given 
threshold uthie as follows: 


= 


ICp(xl^thre) pL 




J r{x)'^ 

where the window function Wp(x) is given by 


Il^p(xl^thre) — 


1 


dM 


dV dn 


iVpeak(^^thre) J d^dfl diVf 

X S{z,M\uthre)bh{z,M), 


(34) 


(35) 


where bh is the linear halo bias. Similarly, the cross power 
spectrum between peaks and convergence is given by (also, see 
Oguri&Takada (2011)) 




PpUe) + Pp!:{£), 


W^peak (^^thre) 

dn 


dz 


d^V 

dadH 


W4x) 


dM S(z, Ml^thre) ^ , > ( \2 

dM pm[z) r[x) 


(36) 


Pp!:{£) = 


dx 


XMUm \k = 


ILp(xkthre)IL4x) 


■(x) 


2(x),M .(37) 


P(X)2 


xfL, ( k jZ{x) 


(38) 


where Mum is the Fourier transform of NEW profile. 

In this paper, we adopt the functional form of bh proposed by 
Bhattacharya et al. (2011) for wCDM model. For HS models, 
we derive the linear halo bias based on the peak-background 
split formalism (Sheth & Tormen 1999). For this purpose, we 
use the model of mass function shown in Li & Hu (2011) and 
relate the linear halo bias with mass function as follows: 


bh{z,M) = 1 


ain/MF 

dSc 


(39) 


where Sc is the linear critical density for spherical collapse 
model, /mf is defined by 




OlnCT, 


ainM 


(40) 


where ctm is the mass variance for HS models. The functional 
form of /mf and om are found in Li & Hu (2011). 


3.4 Dependence of cosmological model 

Figure 1 summarizes the dependence of the weak lensing statis¬ 
tics on cosmological parameters. In order to calculate the lens¬ 
ing statistics for different cosmological models, we use the de¬ 
scription summarized in Section 3.3 and Table 1. The four 
panels in Figure 1 show Ppp, Pp„ and A^peak for vari¬ 
ous cosmological models. In each panel, the black line shows 
the model prediction for ACDM model, whereas the colored 
line represents the model for wCDM models or f{R) mod¬ 
els. In this figure, we adopt the value of wo = —0.8 and —1.2 
or |/flo| = 10“"^ and 10“®. Other cosmological parameters 
are consistent with those derived by WMAP nine-year analy¬ 
sis (Hinshaw et al. 2013). We also present the ratio of each 
statistic with respect to the ACDM model, i.e., we plot 
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multipole I ^thre 

Fig.1. The dependence of weak leasing statistics on model parameters. We plot the power spectra of convergence (top-left), the power spectra of lensing 
peaks (top-right), the cross spectra between convergence and lensing peaks (bottom-left), and the number counts of lensing peaks (bottom-left). We define 
the lensing peaks as those with signal-to-noise ratio r^thre > 3 in the top right and the bottom left panel. The error bars in each panel represent the Gaussian 
or Poission errors in the case of sky coverage of 20,000 square degrees. The dashed line in two top panels corresponds to the shot noise term of auto power 
spectrum. In each panel, we show the theoretical model of weak lensing statistics for cosmological models with different value of wq and |/ho I- The bottom 
portion in each panel provides the ratio of each statistic between ACDM model and «iCDM model or f{R) model. 


(41) 

-AaCDM 

where X = P„„,Ppp,PpK, A^peak. In Figure 1, we simply con¬ 
sider the Gaussian covariance of power spectra and the Poisson 
error of the number count of lensing peaks assuming a hypo¬ 
thetical survey with sky coverage of 20,000 squared degrees. 

The wCDM model mainly affects the linear growth rate 
of matter density perturbations and the cosmic expansion rate. 
Thus, it modifies in an almost scale-independent way and 
causes larger impact on number of more massive halos, i.e. 
IVpeak for higher peaks. 

On the other hand, the f(R) model induces the scale- 
dependence of the linear gravitational growth and reduces grav¬ 
itational force enhancement as a function of the depth of the 
gravitational potential wells of dark matter halos. The scale- 
dependence of growth of linear density fluctuations is deter¬ 


mined by the mass of additional scalar degree of freedom fu in 
this model. Essentially, the fn field characterizes two physical 
length scales: One is the region where the linear growth rate is 
quite similar to that of ACDM and the other is the region where 
the fluctuation growth is amplified by the additional fifth force. 
The transition between the two would occur at the compton 
wavelength of fn. When density fluctuations grow, another im¬ 
portant mechanism, the so-called chameleon mechanism, oper¬ 
ates and induces the density dependence of the mass of fn field. 
Even at the small length scales where the linear growth rate is 
amplified, the chameleon mechanism can effectively switch the 
gravitational force so that ACDM-like evolution is recovered. 
This additional effect would make matter power spectrum and 
halo mass function more complex. Especially, specific features 
like bumps are found in the matter power spectrum and in halo 
mass function on f{R) model (Li & Hu 2011; Zhao 2014). 
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Therefore, the f(R) model amplifies in a wavenumber de¬ 
pendent manner, and changes appreciably the number of halos 
with masses near the relevant mass scale of the force enhance¬ 
ment. 

Although the error estimate in Figure 1 might seem opti¬ 
mistic, Pkk and Apeak are expected to be useful to extract cos¬ 
mological information of the deviation from ACDM model. 
Note, however, that the statistics also depend on parameters 
other than wq and /ro- For example, the convergence power 
spectrum has strong degeneracy among cosmological parame¬ 
ters, while the number of lensing peaks can be much depen¬ 
dent of the concentration of dark matter density profile. In 
order to overcome the degeneracies and uncertainties, we uti¬ 
lize the cross power spectrum of convergence and lensing peaks 
PpK*- The cross correlation can be used to distinguish two dif¬ 
ferent information contents, i.e., cosmological parameters and 
the dark matter density profile. The former can be mainly ex¬ 
tracted from the cross correlation at degree scales whereas the 
latter can be derived from the cross correlation signals at ar- 
cmin scales. Therefore, combined analysis with Pkk, Apeak and 
Pp„ is a powerful tool of cosmology with weak lensing selected 
clusters. 


4 FORECAST 


We perform a Fisher analysis to make forecast for constraints on 
wCDM model or f{R) model with the ongoing Hyper Suprime- 
Cam (HSC) survey. 

Let us briefly summarize the Fisher analysis. For a multivari¬ 
ate Gaussian likelihood, the Fisher matrix Fij can be written as 

Fij = ipr [AiAj + , (42) 

where Ai = C~^dCjdpi, Hij = 2 (dp/dpi) {dp/dpj), C is the 
data covariance matrix, p represents the assumed model, and p 
describes parameters of interest. In the present study, we con¬ 
sider only the second term in Eq. (42). Because C is expected 
to scale proportionally inverse to the survey area, the second 
term will be dominant for a large area survey (see, e.g., Eifler, 
Schneider & Hartlap (2009)). We consider the following set 
of parameters to vary: (10'^As,ns,nnio,Avir,o,tcoor|/Ro|) 

where Avir.o represents the normalization of concentration of 
dark matter density profile. We define Avir,o as 


— Av 


M 


10i4/i-iMo 


(1 + ^)- 


(43) 


'* In principle, the auto power spectrum of lensing peaks Ppp can be a direct 
measure of clustering of clusters. However, we find that the signal to noise 
ratio of Ppp is significantly smaller than that of Pp,^. In the case of lensing 
surveys with the sky coverage of ~1000 squared degrees, we can find 
only ~ 100 peaks at most. Thus, the auto-power spectrum Ppp would be 
dominated by the Poisson noise. We expect that Ppp can be useful to 
extract some information oniy if the survey area exceeds 0,000 squared 
degrees 



10 ® A« 

TLs 

^mO 

44.vir,0 

Wq 

10®1/ro| 

fiducial 

2.41 

0.972 

0.279 

5.72 

-1 

1 

dp 

0.1 

0.13 

0.01 

1.0 

0.2 

0.5 


Table 2. Parameters for Fisher analysis. We denote by dp the 
variation of parameters around the fiducial vaule that is used for 
estimation of the numerical derivative of Pkk, PpK and Apeak- 


We construct the data vector D from a set of binned spectra 
Pkk,Ppk and the number count of peaks Apeak as, 

Di = {P^Kiil) , ■■■, PKK{ilo) , Ppniil), ■■■, Ppniiio), 

Apeak(Pthre,l), Apeak(Pthre,6)}7 (44) 

where Ppa is defined by the cross spectrum between conver¬ 
gence and lensing peaks with Uthre = 3. For the Fisher analysis, 
we use 10 bins in the range of £i = [20, 2000] and 6 bins in 
the range of Pthre,i = [2.5,5.5]. In total, a data vector has 26 
elements, 2 x 10 for power- and cross-spectra and 6 for peak 
counts. In the Fisher analysis, we need the derivative of statis¬ 
tics Di with respect to p. We evaluate the numerical derivatives 
as follows: 

dPj _ Dj {p^a^ +dpa)- Dj (pi°^ - dp a) 

dpa 2dpa 

where p^a^ is the fiducial value, dpa is the variation of a-th pa¬ 
rameter and we use the similar definition to Pp^ and Apeak- 
In order to evaluate the term dDi jdpa in Eq. (45), we calcu¬ 
late the relevant statistics adopting the analytic models shown 
in Section 3.3 and Table 1 for different cosmologies. We sum¬ 
marize the fiducial value of p and dp in Table 2. 

We need a 26 x 26 data covariance matrix for the Fisher anal¬ 
ysis. In order to derive the covariances, we use 200 masked sky 
simulations genearated in Paper I. 

The masked sky simulations are constructed based on ten in¬ 
dependent full-sky weak lensing simulations. We perform the 
full-sky weak lensing simulations combined with a large set of 
dark matter A-body simulations with consistent cosmological 
model with WMAP nine-year results (Hinshaw et al. 2013). In 
each simulation, we properly take into account the deflection of 
light path on a curved sky. From each full-sky realization, we 
extract multiple “observed” regions with similar survey geom¬ 
etry to the Hyper-Suprime Cam (HSC) survey. We create 20 
convergence maps on masked HSC sky by selecting the desired 
sky coverage (565 -i- 680 = 1245 squared degrees) from single 
full-sky simulation. We incorporate the masked regions associ¬ 
ated with the positions of bright stars inside the HSC survey re¬ 
gions. The details of the full-sky weak lensing simulations and 
masked regions are found in Paper I. From 200 realizations of 
mock HSC data, we construct smoothed convergence maps with 
the compensated Gaussian filter and fully utilize these maps to 
derive the covariance matrices between the statistics. 

When calculating the inverse covariance, we include a de- 
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biasing correction as a = (rireai — n-bin — 2 )/(nreai — 1 ) with 
rirea = 200 being the number of realization of simulation sets 
and ribin = 26 being the number of total bins in our data vector 
(Hartlap, Simon & Schneider 2007). 

We also take into account the constraints from the CMB 
priors expected from the Planck satellite mission. When we 
compute the Fisher matrix for the CMB, we use the Markov- 
Chain Monte-Carlo (MCMC) engine for exploring cosmologi¬ 
cal parameter space COSMOMC (Lewis & Bridle 2002). We con¬ 
sider the parameter constraints from the angular power spec¬ 
tra of temperature anisotropies, iJ-mode polarization and their 
cross-correlation. For the MCMC analysis, in addition to 
Iff As,ns,Q.ch? and wq, our independent variables include the 
baryon density Q.bh^, Flubble parameter h, reionization optical 
depth T, To examine the pure power of weak lensing statistics 
to constrain wo and |/ro|, we do not include any constraints 
on Wo and |/ro| from the CMB. Assuming that the constraints 
from the CMB and the lensing statistics are independent, we 
express the total Fisher matrix as 

F = Flensing + FcMB . (46) 

When we include the CMB priors by Eq. (46), we marginal¬ 
ize over the other cosmological parameters except the following 
three parameters; if Aa,ns, and flmo- 

Figure 2 summarizes the expected constraints on two cos¬ 
mological models (wCDM or f{R)) in the case of HSC survey. 
The left panel shows the marginalized Icr constraints on wCDM 
models in two-dimensional parameter space and the right panel 
shows the constraints on f{R) models. In each panel, we in¬ 
dicate constraints (error circles) from different statistics by col¬ 
ored lines. 

In the case of wCDM model, there remains strong degen¬ 
eracy among parameters; especially, it is difficult to break the 
degeneracy among wo. As and Dmo with our statistics alone. 
This is because we can measure the amplitude of matter density 
fluctuations by using both Pk,k, and Apeak and there apprear no 
characteristic features in our lensing statistics. In order to break 
such degeneracies, we require other information such as CMB 
measurements. When including the information from CMB, we 
can successfully break the degeneracy among cosmological pa¬ 
rameters and thus measure wo with a level of ~ 0 . 1 . 

The situation becomes slightly different in the case of f{R) 
model. The modification of gravitational force due to the addi¬ 
tional scalar degree of freedom induces new characteristic fea¬ 
tures, which can not be compensated by other cosmological pa¬ 
rameters. These features are expected to show up as excess of 
the halo mass function at M ~ 10^^ h~^MQ (Li & Hu 2011), or 
excess in the lensing power spectrum at multipole £ ~ 100 — 500 
(Zhao 2014). The relevant mass and length scales of the ex¬ 
cess are set by the value of |/ro|. Our weak lensing statistics 
presented here are sensitive and thus one can measure such in¬ 


teresting features using data from upcoming lensing surveys. 
Interestingly, the expected constraint on |/ro| is at a level of 
10“® with our weak lensing statistics alone. In addition to our 
lensing statistics, parameter constraints from CMB measure¬ 
ments would greatly help to break degeneracies among other 
parameters than |/ro| and Avir.o. Since the CMB information 
provides generally more stringent constraints on Us and As, we 
can further improve the constraint on | /ro | by combining CMB 
with our lensing statistics. The expected Icr error for each pa¬ 
rameter is summarized in Table 3. 

5 CONCLUSION AND DISCUSSION 

We have studied cosmological information content in statistics 
of weak-lensing selected clusters. We have developed a theo¬ 
retical model of the lensing statistics and applied to two com¬ 
peting cosmological models: wCDM model and f{R) model. 
The key parameter in wCDM model and f{R) model are the 
equation of state parameter of dark energy wo and the addi¬ 
tional scalar degree of freedom |/ro|, respectively. We have 
utilized masked sky simulations (see Paper I for details) to es¬ 
timate non-Gaussian covariance caused by non-linear gravita¬ 
tional growth and the mode-coupling due to masked regions 
simultaneously. With such non-Gaussian covariance, we have 
performed a Fisher analysis, which yields realistic forecast for 
constraining the nature of dark energy or the modification of 
gravity. 

We consider specifically ongoing Hyper Suprime-Cam 
(HSC) survey with a sky coverage of ~ 1250 squared degrees. 
Combined analysis of cosmic shear and weak lensing selected 
clusters can constrain wo with a level of Awo ~ 0.1, and |/ro| 
with a level of A|/ro| ~ 5 x 10 “® with the help of cosmic 
microwave background measurements. Note that the expected 
constraint on | /ro| is comparable to the recent constraints of 
(Cataneo et al. 2014). Clearly, our approach is promising for 
upcoming surveys with a sky coverage of 10,000 squared de¬ 
grees or more. Assuming that the statistical error in upcom¬ 
ing wide field surveys is reduced proportionally to the effec¬ 
tive survey area, we can improve the constraints by a factor 
of ~ (20000/1250)^/^ = 4 in the case of the Large Synoptic 
Survey Telescope (LSST)^ with a proposed sky coverage of 
20,000 square degrees. We summarize the expected error with 
a LSST-like survey in Table 4. 

The normalization of power spectrum of matter density per¬ 
turbations is one of the basic quantities in measure of gravi¬ 
tational growth at low redshift. In practice, we normalize the 
linear matter power spectrum by using the following quantities: 

/ p j3 \l/2 

aRiz)=( -^\WTH{kR)\^P^ik,z)\ , (47) 

^ http://www.lsst.org/lsst/ 
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wCDM 

10 ® A, 

ris 

^mO 

4lvir,0 

Wo 

(HSC) 

1.43 (0.0409) 

0.111 (0.00385) 

0.0918 (0.00192) 

-(-) 

0.686 (0.0227) 

P.. + Pp. (HSC) 

1.40 (0.0406) 

0.109 (0.00381) 

0.0897 (0.00190) 

1.77 (1.58) 

0.676 (0.0225) 

+ A^peak (HSC) 

0.548 (0.0338) 

0.0830 (0.00331) 

0.0728 (0.00170) 

1.14(0.314) 

0.469 (0.0193) 

Pkk + P-pK + A^peak (HSC) 

0.497 (0.0299) 

0.0758 (0.00291) 

0.0669 (0.00146) 

1.02 (0.306) 

0.417 (0.0165) 

Pkk + PpK + A^peak (HSC) + Plaiick 

0.132 (0.0294) 

0.0145 (0.00274) 

0.0124 (0.00143) 

0.379 (0.306) 

0.108 (0.0165) 

f{R) 

10 ® A, 

ris 

^mO 

^vir,0 

10®|/ro| 

P.. (HSC) 

1.35 (0.0386) 

0.268 (0.00774) 

0.0951 (0.00170) 

-(-) 

5.91 (0.180) 

+ (HSC) 

1.29 (0.0384) 

0.224 (0.00771) 

0.0768 (0.00169) 

1.82 (1.65) 

4.79 (0.178) 

Pk.k. + A^peak (HSC) 

0.571 (0.0309) 

0.0168 (0.00372) 

0.0311 (0.00146) 

1.46 (0.296) 

0.897 (0.129) 

Pk.k. + PpK + A^peak (HSC) 

0.418 (0.0267) 

0.0136 (0.00297) 

0.0219 (0.00133) 

0.913 (0.291) 

0.717 (0.127) 

Pkk + PpK + A^peak (HSC) + Plaiick 

0.0898 (0.0263) 

0.00848 (0.00278) 

0.00735 (0.00130) 

0.677 (0.291) 

0.481 (0.127) 


Table 3. The summary of the expected marginalized error with la confidence level. In order to make the forecast of constraints 
in this table, we assume the following parameters: the rms of intrinsic ellipticies of sources a^y = 0.4, number density of sources 
rigai = I0arcmin“2^ and source redshift Zaource = 1. For comparison, we also provide the unmarginalized error as the number in brackets. 


uiCDM 

10 ® A, 

ris 

^mO 

-^vir,0 

Wo 

(LSST) 

0.359 

(0.0102) 

0.0279 

(9.64xl0-'‘) 

0.0229 

(4.80xl0-'‘) 

(-) 

0.171 

(0.00568) 

Pkk + PpK (LSST) 

0.350 

(0.0101) 

0.0272 

(9.53x10"^) 

0.0224 

(4.76x10"^) 

0.443 

(0.395) 

0.169 

(0.00564) 

Pkk. + A^peak (LSST) 

0.137 

(0.00845) 

0.0207 

(8.29x10-^) 

0.0182 

(4.26x10-^) 

0.286 

(0.0785) 

0.115 

(0.00483) 

Pkk. + PpK “h A^peak (LSST) 

0.124 

(0.00748) 

0.0189 

(7.28xl0"'‘) 

0.0167 

(3.67xl0-'‘) 

0.256 

(0.0767) 

0.104 

(0.00414) 

Pkk + PpK. + A^peak (LSST) + Plaiick 

0.0818 

(0.00747) 

0.0115 

(7.25x10-^) 

0.0101 

(3.66x10-^) 

0.169 

(0.0767) 

0.0665 

(0.00414) 

fip) 

10 ® A, 

ris 

^mO 

-^vir,0 

10®|/ro| 

(LSST) 

0.339 

(0.00965) 

0.0671 

(19.3xl0-'‘) 

0.0237 

(4.26xl0-'‘) 

(-) 

1.47 

(0.0452) 

Pkk + PpK (LSST) 

0.322 

(0.00960) 

0.0562 

(19.2xl0“'‘) 

0.0192 

(4.22xl0“'‘) 

0.455 

(0.414) 

1.19 

(0.0447) 

Pkk + A^peak (LSST) 

0.142 

(0.00773) 

0.00422 

(9.30xl0-'‘) 

0.00779 

(3.65x10-^) 

0.366 

(0.0741) 

0.224 

(0.0324) 

Pkk + PpK + A^peak (LSST) 

0.104 

(0.00668) 

0.00340 

(7.41xl0"'‘) 

0.00586 

(3.32x10-^) 

0.228 

(0.0729) 

0.179 

(0.0319) 

Pkk + PpK + A^peak (LSST) + PlailCk 

0.0655 

(0.00667) 

0.00285 

(7.38xl0-'‘) 

0.00369 

(3.32xl0-'‘) 

0.204 

(0.0729) 

0.151 

(0.0319) 


Table 4. The expected marginalized la error in a LSST-like survey with sky coverage of 20,000 square degrees. 
We adopt the same parameter as shown in Table 3. We scale the covariance matrices with survey area (i.e. 
by a factor of 20000/1250). As for Table 3, the unmarginalized error is shown as the number in brackets. 
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with non Gaussian error _2 — 1.5 —1 —0.5 0 with non—Gaussian error 0 5 10 

Fig. 2. The error circle obtained from Fisher analysis in the case of Hyper Suprime-Cam survey. In this figure, we show the cosmological constraints from 
weak lensing statistics alone. Each colored line corresponds to the constraints derived from different statistics: (red), and Pp„ (green), P„k and 

fVpeak (blue), and combined with three (cyan). The left panel represents the constraint on u;CDM model, while the right one shows the constraint on f(R) 
model. To estimate the constraints presented here, we take into account the non-Gaussian covariances derived from 200 masked sky simulations. 


where WTH{kR) is the Fourier transform of top hat func¬ 
tion with scale of R and R is commonly set to be 8 Mpc/h. 
Theoretically, this quantity 0-8(2) can be derived when the lin¬ 
ear matter power amplitude is measured at higher redshift and 
the linear growth rate is taken into account properly. On the 
other hand, we can measure the value of 0-8(2) in a more di¬ 
rect manner when we probe the clustering of matter density 
field at lower redshift. Therefore, comparing the value of as 
derived from high redshift information with the observed as at 
lower redshift provides an invaluable opportunity to check con¬ 
sistency of cosmological models between high and low redshift 
(MacCrann et al. 2014). Figure 3 shows the expected constraint 
on as with the lensing statistics proposed in this paper. In this 
figure, we present the expected constrained region on flmo — 
as{z = 0) plane with 95% confidence level. In order to derive 
such regions, we sample 100,000 points in five cosmological 
parameters space for uiCDM model (i.e. As, ris, flmo, 4lvir,o 
and wq) assuming Gaussian distribution function with covari¬ 
ance estimated from Fisher analysis (see. Section 4 for details). 
Then, we calculate the value of 0-8(2 = 0) at each point. With 
prior information of CMB measurement, our statistics can con¬ 
strain on 0-8(2 = 0) with a level of 0.025 and 0.01 for HSC-like 
and LSST-like survey, respectively. Interestingly, in the case of 
LSST-like survey, we find that 0-8(2 = 0) can be constrained 
with precision of ~0.025 even in absence of CMB informa¬ 
tion. In Figure 3, the square and star symbols represent the de¬ 
rived values of 0-8 for wCDM and f{R) models for a given As. 
For As derived from CMB, dark energy and modified gravity 
model predict the different values of 0-8 from that of the stan¬ 
dard ACDM model. Thus, the constraint on 0-8 by our lensing 
statistics serves as an important test on the self-consistency of 


0.95 


0.9 


0.85 


0.8 



0.26 0.28 0.3 0.32 


Fig. 3. The expected constraints on fimo - as plane. The green and white 
region correspond to the marginalized 95% confidence level for HSC and 
LSST-like survey, respectively. We evaluate these regions from the Fisher 
matrix in combined analysis of cosmic shear power spectrum, the cross 
correlation of convergence and its peak, and the number count of peaks 
with prior information of CMB measurement. The center of confidence 
surface represent the standard ACDM model, while square and star 
symbols show the different value of o-g in ujCDM model and f{R) model. 


the standard ACDM model through 2 ~ 1000 to 2 ~ 0. 

Combined statistical analysis with weak lensing selected 
clusters and cosmic shear provides precise measurement of 
gravitational growth at low redshift, and thus can place robust 
constraints on variant cosmological models. However, in order 
to apply our method to real data set, we need to consider sys¬ 
tematic effects in detail. One of the most important systematics 
is source redshift uncertainty. Throughout this paper, we have 
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assumed all the source galaxies are located at the same redshift. 
In reality, source galaxies are distributed over a wide redshift 
range. In order to examine the effect of source redshift uncer¬ 
tainty, we compare the case of ^source = 0.9 and ^source = 1 un¬ 
der ACDM cosmology. The 10% difference in the mean source 
redshift roughly causes under/over-estimate of flmo with a level 
of ~ 0.01 for cosmic shear power spectrum and cross power 
spectrum of convergence and its peaks, while the number count 
of lensing peaks is less sensitive. Therefore, the mean source 
redshift needs to be calibrated with a level of 0.1 for HSC sur¬ 
vey. Clearly, it is important to study the effect of photometric 
redshift uncertainty in a more quantitative manner, using, for 
instance, realistic mock galaxy catalogs. 

The baryonic processes such as gas cooling and stellar feed¬ 
back cause bias in parameter estimation with our lensing statis¬ 
tics. Previous studies (Semboloni et al. 2011; Semboloni, 
Hoekstra & Schaye 2013; Zentner et al. 2013) indicate that 
the baryonic effect on matter clustering could change two-point 
statistics of cosmic shear at £ ~ 1000 by a few percent. Since 
a 3% difference of roughly corresponds to Aflmo ~ 0.015, 
the baryonic effect needs to be accounted for in the case of 
LSST-like survey. For peak statistics, we already take into 
account a sort of the baryonic effect by considering the un¬ 
certainty of halo concentration. We argue that the baryonic 
processes would cause minor effect on peak statistics for high 
peaks, as Osato, Shirasaki & Yoshida (2015) have examined. 
Nevertheless, the detailed modeling of baryonic effect on weak 
lensing selected clusters will allow to derive more robust con¬ 
straint on cosmological models. 

Massive neutrinos are expected to have a significant effect 
on cosmic structure formation (e.g.. Bond, Efstathiou, & Silk 
1980). According to Harnois-Deraps et al. (2015), massive 
neutrinos with Mi, — 0.2 eV suppress the amplitude of con¬ 
vergence power spectrum at ^ ~ 1000 by a factor of ~ 0.9 
(see Figure 1 and 2 in Flarnois-Deraps et al. (2015)). In ad¬ 
dition, Costanzi et al. (2013) have shown massive neutrinos 
with Ml, = 0.2 eV change the halo mass function with the halo 
mass of ~ lO’^'^ h~^MQ by a factor of ~ 0.7. Therefore, ignor¬ 
ing massive neutrinos can potentially cause the mis-estimation 
of (78 and/or flmo- Further studies on the impact of massive 
neutrinos on weak lensing statistics are warranted. 

Magnification by weak lensing causes the scatter of the 
brightness and/or the size of source galaxies. This magnifica¬ 
tion effect on cosmic shear statistics has been investigated in 
literature (Schmidt et al. 2009b; Schmidt & Rozo 2011; Liu 
et al. 2014). For convergence power spectrum, Schmidt et al. 
(2009b) have estimated that the magnification effect typically 
induces ~1% differences of at 1000. On the other hand, 
the peak height by the NFW dark matter halo is less sensitive 
to the magnification effect even if one reconstructs the lensing 
mass map with a smoothing scale of a few arcminutes (e.g., see 


Figure 5 in Schmidt & Rozo (2011)). The magnification effect 
on peak statistics would be important in the case of LSST-like 
survey as Liu et al. (2014) pointed out. 

Other than magnification effect, there exist possible system¬ 
atic effects which could be critical in cosmological analysis with 
weak lensing selected clusters. Source-lens clustering (Flamana 
et al. 2002) and the intrinsic alignment (Flirata 8l Seljak 2004) 
are likely to compromise cosmological parameter estimation. 
Although the impact of those effects on cosmic shear statistics 
is still uncertain, a promising approach in theoretical studies 
would be associating the source positions with their host dark 
matter halos on the light cone. 

Future galaxy imaging surveys will generate a large amount 
of high-precision cosmic shear data. It will then become possi¬ 
ble to map the matter density distribution in the universe and to 
measure the growth of the large-scale structure in a direct man¬ 
ner. The lensing statistics studied in this paper will enable us 
to make the best use of data from future large surveys to ex¬ 
plore the origin of cosmic acceleration, or the deviation from 
the concordance ACDM cosmology. 
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